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Abstract: We study lepton flavor- violating (LFV) processes within a supersymmetric 
type-I seesaw framework with flavor-blind universal boundary conditions, properly ac- 
counting for the effect of the neutrino sector on the dark matter relic abundance. We 
consider several possibilities for the neutrino Yukawa coupling matrix and show that in 
regions of SUSY parameter space that yield the correct neutralino relic density, LFV rates 
can differ from naive estimates by up to two orders of magnitude. Contrary to common 
belief, we flnd that current LFV limits do not exclude neutrino Yukawa couplings larger 
than top Yukawa couplings. We introduce the ISAJET-M program that was used for the 
computations. 



1. Introduction 



It is now a firmly established experimental fact that neutrinos are massive and mix Q. A 
global fit to current neutrino oscillation data |^] gives the following 3a ranges/limits for 
the mixing parameters (in the "standard parametrization" , Eq. A.14| ) , 
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The phases are not constrained by current data. 

The most elegant and popular explanation for small neutrino masses is offered by the 
seesaw mechanism Q. Here the Standard Model (SM) is augmented by three right-handed 
neutrinos (RHNs). The RHNs transform as singlets under the SM gauge group and thus 
can have a Majorana mass, which is taken to be very large, Mpjaj ^ GeV. One of 
the implications of neutrino mixing is the possibility for charged lepton flavor-violating 
processes. In the SM-seesaw these rates are highly suppressed by the very large scale 
MMaj Hi) so that LFV processes are a good probe of new physics. 

Supersymmetry (SUSY) is a well motivated possibility for new physics [^, ^ because, 
among other things, it stabilizes the seesaw mechanism [0] and ameliorates the hierarchy 
problem appearing from the presence of the very high scale Mm 
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SUSY must be 

a broken symmetry that is parametrized by soft SUSY breaking (SSB) terms in the La- 
grangian . In general these terms can have arbitrary flavor structures that would induce 
unacceptably large flavor violating effects - this is the well-known SUSY flavor problem. 
The simplest solution is to assume a flavor-blind SUSY breaking mediation mechanism 
that will generate flavor universal SSB terms at some high scale. However, this does not 
mean that SSB terms will also remain flavor-universal at the weak scale: flavor-violating 
terms will be generated by the Yukawa couplings during the evolution from the high to the 
weak scale. As a consequence, SUSY contributions to LFV processes are suppressed by the 
characteristic mass scale of SUSY particles MsusY ~ 1 TeV (instead of Mpjaj) and thus 
LFV can be observable. Many authors have studied these processes under various SUSY 
model assumptions and seesaw parameters (see e.g., |lO|, |ll|, ||, ||, |l5|, 0, ||, ||]). 

The prediction of LFV rates requires knowledge of the neutrino Yukawa coupling ma- 
trix. However, experimentally measured light neutrino masses and mixings do not provide 
sufficient information to determine it [0]. A top-down approach is frequently adopted in 
which the neutrino Yukawa matrix is set by a specific SUSY-Grand Unified Theory (GUT), 
with an 50(10) gauge group being the favorite choice. In this work we take a general ap- 
proach, and consider two cases of Yukawa unification parameters (defined in Eq. 2.13| ) 
inspired by 5*0(10) relations. For both of these cases we study two scenarios with small 
and large mixings in the neutrino Yukawa matrix. Our analysis is based on the type-I 
seesaw mechanism Q. 

The existence of a massive, electrically and color neutral, stable weakly interacting par- 
ticle that can serve as a cold Dark Matter (DM) candidate is perhaps the most compelling 
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feature of the i?— parity conserving Minimal Supersymmetric Standard Model (MSSM). In 
most cases the DM particle is the lightest neutralino, Z\ |2^. The mass density of DM has 
been precisely determined by cosmological measurements: a combination of WMAP CMB 



data with the baryon acoustic oscillations in galaxy power spectra gives [21 



where = pj Pc with pc the critical mass density of the Universe, and h is the scaled Hubble 
parameter. Such a precise determination places severe constraints on new physics scenarios. 
In the simplest SUSY model with universal SSB values at the high scale, mSUGRA (or 
CMSSM) ||2^, only a few regions of parameter space survive: the bulk region |2^, the 
stau |p5| , or stop [27| coannihilation regions, the hyperbolic branch/focus point (HB/FP) 
region |28|, and the A ox h resonance annihilation (Higgs funnel) regions |24, ^9|, |3^. The 
benchmark values of mSUGRA input parameters for these regions are listed in Table [l|. 

The introduction of the RHNs 
and their associated Yukawa cou- 
plings changes predictions for some 
sparticle masses via new contribu- 
tions to the renormalization group 
equations (RGEs). These imprints 



couplings in collider experiments |31]. 
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Table 1: Input parameters for benchmark points and 
corresponding DM-allowed regions of mSUGRA. The di- 
mensionful parameters mg, ?ni/2 and are in GeV. For 
all points > and = 171 GeV. 



In previous work [32|, some of us 
demonstrated that these neutrino- 
induced changes in the SUSY mass 
spectrum can significantly alter the 
DM ( CO ) annihilation mechanisms with 
concomitant changes in direct and indirect detection rates. 

The aim of the present work is to study predictions for LFV rates while correctly taking 
into account the aforementioned effect on the neutralino relic density. We take a model- 
independent approach and only consider effects from RGE running below the unification 
scale Mgut- These two important points distinguish our work (which is closer in spirit to 
the study of Ref. 0) from the study of Ref. [||. We find that proper consideration of 



the interplay between the neutrino and SUSY sectors can change the predictions for the 
LFV rates in WMAP-allowed regions by a factor between 2 — 100 compared with naive 
estimates, and that contrary to common belief, large neutrino Yukawa couplings are not 
ruled out by current LFV bounds. 

The rest of the paper is organized as follows. In the next section we briefly review 
LFV processes in SUSY-seesaw framework and motivate our ansatz for the neutrino Yukawa 
matrix. Results from numerical analyses are presented in Section ^. Section ^ is devoted 
to the discussion of our findings. Finally, we present our conclusions in Section |5[ Our 
notation, a description of our code, and Yukawa RGEs with thresholds can be found in the 
appendices. 
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2. SUSY-seesaw and LFV processes 



We begin with a brief discussion of our formalism. Details about our notation and con- 
ventions are relegated to Appendix The superpotential for the MSSM augmented by 
singlet right-handed neutrinos Nf is 



/ 



fhlSSM + (fj/) 



(2.1) 



where Jmssm is the MSSM superpotential shown in Eq. ( |A.1| ), L and Hu are, respectively, 
lepton doublet and up-higgs superfields, and M^v represents the Majorana mass matrix for 
the (heavy) right-handed neutrinos. Above the scale Muaj the light neutrino mass matrix 
is given by the well-known type-I seesaw formula 



(2.2) 



where is the vacuum expectation value (VEV) of the neutral component /i^ of the up- 
type Higgs doublet H^- At low energies the RHNs decouple from the theory and the light 
neutrino mass matrix is A^j^ = — ku^. Here k is a coupling matrix of the dimension-5 effec- 
tive operator generated by RHNs (see Eq. ( A.ll| ) for the definition), which is determined 
by matching conditions ( A.12| ) at the RHN decoupling thresholds. 
A common solution to the SUSY 



flavor problem is to assume that 
at some high scale sfermion SSB 
mass-squared matrices are diagonal 
and universal in fiavor, and the tri- 
linear couplings are proportional to 
the Yukawa couplings. In the frame- 
work of mSUGRA extended by right- 
handed neutrinos (mSUGRA-seesaw) 
the SSB boundary conditions at Grand- 
Unification Scale Mgut take the Table 2 
particularly simple form: 
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Since the Yukawa couplings fg and fi, cannot be simultaneously diagonalized, non-vanishing 
off-diagonal elements will be generated in the left-handed slepton mass matrix mj^ via 
renormalization group evolution. In the leading-logarithmic approximation with universal 
boundary conditions (^.3]) the off-diagonal elements are 
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(mi),^, -^(3m2 + Al) ^(fj),fe(f:),, log 



(2.4) 



where M]\j^ are RHN decoupling scales that are approximately equal to the eigenvalues 
of the Majorana mass matrix M^r. At low energies, these off-diagonal terms induce LFV 



processes such as /, — > /j7, U 



3lj and /j 



Ij conversion in nuclei. 



Current bounds on 



LFV processes as well as the projected sensitivities of future experiments are summarized 
in Table |. 
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The branching ratio for the flavor- violating radiative decay of a charged lepton is given 

by 

BRik - /,-7) = + \Ar\') ■ (2.5) 

Here a is the electromagnetic fine structure constant, T{li) is the total decay width of the 
initial lepton, and A^ ji are form factors for left and right chiralities of the incoming lepton 
whose full expressions in SUSY were obtained in Ref. il^. Because m;. » m,;. one has 



Ar ^ Al in the case of initial universality such as (|2.3| ) llQ , 13 1. In the mass-insertion 
approximation, the branching ratio can be related to the corresponding off-diagonal element 
of the left-handed slepton mass matrix [O], 



a3 



BR{li ^ Ij-f) ~ BR{li ^ Iji^ji^i)-^;^ tan^ f3 , (2.6) 

where Gp is the Fermi coupling constant and is the characteristic mass scale of the SUSY 
particles in the loop. In the case of universal boundary conditions (|2.3| ), this expression 
used in conjunction with the leading- log result (vA) well approximates the full expression 



BJ), if one sets [14| 



LFV li — > 3lj decays and /j — > L conversion occur via 7-, Z— and Higgs-penguins as 



well as squark/slepton box diagrams |12]. Higgs-penguins dominate in the regime of large 
tan/3 ~ 60 and light H/A Higgs boson mass ~ 100 GeV, and enhance rates by up to a few 
orders of magnitude [|l^]. However, the latter condition cannot be realized in the universal 



scenario (2.3) |4^]. It was shown in Ref. [^] that these LFV rates are well described in the 
universal scanario by the same 7-penguins that contribute to the radiative decays. The 
branching ratio for trilepton decays is approximately given by 

BR{k ^ 31,) ^ ^ (log ^ - ^ ) BR{k ^ Z,-7) • (2.8) 
37r \^ mf^ 4 J 

For ^ e conversion a similar relation holds: 

CRifiN ^ eN) ^ ^ = Z^jf \F{q')\' BRit, ^ 67) , (2.9) 

-L capt J- capt 

where Z is the proton number of the nucleus A^, Z(,ff is an effective atomic charge obtained 
by averaging the muon wave function over the nuclear density 1^^, F{q^) denotes the 
nuclear form factor at momentum transfer q [^] and T^apt is the measured total muon 



capture rate |48]. In this work we consider two target materials that will be used by 
future experiments. For 22^1, which will be used by the proposed PRIME experiment at 
J-PARC II, ^e// = 17.6, F(g2 ~ -ml) ~ 0.54 and T^apt = 2.590 x lO'^sec-i. For f^Al, 



the target material for the proposed Mu2e experiment at Fermilab |44|, Z^ff = 11.48, 
F{q^ ~ -ml) ~ 0.64 and T^apt = 7.054 x lO^sec'^. 
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2.1 50(10) GUTs 

As discussed in the previous section, LFV rates crucially depend on the neutrino Yukawa 
coupling matrix iiy. However this matrix cannot be reconstructed from experimental data 
by inverting the seesaw formula ( |2.2D : f^, and M^v together depend on 18 parameters, while 
Aiu contains only 9 observables. A common solution is to turn to GUTs where fj^ is related 
to the known Yukawa matrices of SM fermions. 

SO (10) GUTs unify all SM fermions and the right-handed neutrino of each generation 
in a single 16-dimensional spinor representation. The Higgs representation assignments 
are determined by the following decompositions of the direct products: 

16 16 = 10 ©120 ©126, (2.10) 
160 16 = 1 ©45 ©210. (2.11) 

Many 5*0(10) models exist in the literature with different choices of Higgs represen- 
tations and, frequently, with additional flavor symmetries. These models can be divided 
in two general classes |Q. The first uses only low dimensional Higgs multiplets 10, 16, 
45 with some nonrenormalizable operators constructed from them. This necessarily leads 
to large iZ— parity violation so that these models cannot provide a DM candidate. Mod- 
els in the other class involve 10, 120 or 126 Higgs representations, have renormalizable 
couplings, preserve iZ— parity, and are often referred to as minimal Higgs models. The 
resulting set of sum-rules for the mass matrices are 
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where i-jz are 50(10) Yukawa coupling matrices, are VEVs of the various SU{2)l 
doublets (with Higgs fields residing in 7^ = 10, 126, 120), and vl and Vr are the B — L 
breaking VEVs of the SU{2)l triplet and singlet respectively. In the type-I seesaw, which 
we consider in this paper, vl = and SUSY prevents it reappearance via loop diagrams Q. 

From Eq. ( |2.12| ), if Higgs superfields reside in 10, as they do in the simplest scenarios, 
then the neutrino Yukawa matrix will be identical to the up-quark Yukawa at Mgut- If 
the higgses are in 126, then f^, = — Sf^. A dominant 120 Higgs would lead to at least a pair 
of degenerate heavy up-quarks and thus is phenomenologically excluded. Motivated 
by the above, we introduce a neutrino Yukawa unification parameter as 

iuOcR^u^u, (2-13) 

and consider two cases, ^ \Ruu\ = 1 and \Ruu\ = 3. Note that f^, and f„ need not be aligned: 
the subdominant contributions from other higgs multiplets and/or flavor group structure 

""^ Henceforth, we denote \R,^u\ by R,^u- 
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can lead to different diagonalization matrices. To keep our discussion as simple as possible, 
we consider two extreme cases of the mixing present in f^. 

Large mixing: The measured values of the neutrino mixing angles ( |1.1| ) are con- 
sistent with the so-called tri-bimaximal pattern |^l], ^ , where sin^ 9i2 = ^ , sin^ 623 = 
^, sin^ = 0. Thus it is reasonable to postulate that mixing in the neutrino Yukawa 
matrix at the GUT scale also has a tri-bimaximal form. In other words, we assume that 
the observed MNS mixing matrix arises only from the left-handed rotation matrix (i.e., we 
set ^i/^ = We take a neutrino Yukawa matrix of the form^ 



f, 



fdiag 



(2.14) 
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where s^^ = sin^, = cos^, -s^ = sine/), = cos</> and the parameters X: 



(2.15) 



G [0,2^]. 



This is the simplest generalization of a tri-bimaximal mixing (x = = 0) that allowes CP 
violation |^2|. A fit to experimental values (^]^) reveals that the Harrison-Scott parameters 
X and (j) are restricted to the vicinity of x + </> — nir, with the upper bound on ^13 imposing 
the strongest constraint. Ignoring RGE effects, we invert the seesaw formula and obtain 
the approximate form for the Majorana mass matrix 



Mat ~ diag 



m 



m 



f2 



m 



X Rt 



(2.16) 



We will consider this case as representative of the large mixing scenario. 

Small mixing: For the small mixing scenario, we take the neutrino and up-quark 
Yukawa matrices to be exactly aligned with each other at the GUT scale, 



R-UU^U ) 



(2.17) 



so that neutrino mixing is given by the CKM matrix. Then, in the absence of significant 
RGE magnification effects, the Majorana mass matrix cannot be diagonal. Assuming tri- 
bimaximal mixing in and neglecting the small mixing in fjy we can estimate eigenvalues 
M jv for the normal hierarchy of light neutrinos (niui ^ "^1^2 ^ "^i^a ) 1 
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For the inverted mass hierarchy {m^-^ ~ ^ ^1/3)1 a similar procedure yields 
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^This corresponds to trivial misalignment matrix R = 1 in the Casas-Ibarra parametrization llOl 
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Notice that the largest RHN mass is controlled by the smallest light neutrino mass. 

From Eqs. ( 2.16| )-( p.l9 ), we see that RHNs have a very strong mass hierarchy "quar- 
datic" to the one in up-quark sector: Mat^ : M^^ : M^^ ~ : : m^. For this reason, 
only the spectrum with normal hierarchy of light neutrinos is feasible. A quasi-degenerate 
spectrum {'my-^ ~ rriy^ ~ "^u-i) would require the lightest Majorana mass to be in the 
10^ — 10^ GeV range with significant L-R mixing in the sneutrino sector which is in confict 
with our approximations (see the discussion pertaining to Eq. |A.15 ). Moreover, such light 
Majorana masses make successful thermal leptogenesis impossible. The inverse hierarchical 
case would require the heaviest Majorana mass to be of order 10^ GeV, which is well above 
the GUT scale. This type of spectrum also suffers instabilities under very small changes 
to Mat and RGE evolution [53|. In the rest of the paper we concentrate on the normal 
hierarchy case, which is also favored by GUT model building 



3. Procedure and results 

We extensively modified ISAJET ||5^ by including the neutrino sector and by implementing 
RGE evolution in matrix form to incorporate flavor effects in both the quark and lepton 
sectors. The resulting program ISAJET-M performs RGE evolution in the MSSM aug- 
mented with RHNs at 2-loop level taking into account various thresholds and computes 
sparticle spectra including radiative corrections. The computation of the neutralino relic 
density is done using the IsaReD code and LEV rates are computed using the full one- 
loop formulae from Ref. |12|. A graphical outline of our procedure is presented in Fig. |l|, 



and details of the program are provided in Appendix Q 

In the neutrino sector we employ a "top-down" approach in which and M ^ are inputs 
at the scale Mqut- Physical neutrino masses and mixings are derived results which we 
require to be within the experimental bounds ( |1.1| ). We consider the two cases for neutrino 
Yukawa unification parameter, R^u = 1 and Ryu = 3 that were introduced earlier. For 
each case we consider scenarios with large and small neutrino mixings using Eqs. ( 2.14| ) and 



( 2.17 ) to set the neutrino Yukawa matrix at Mqut- This also restricts Majorana masses 
to be below that scale, i.e. max(M7vJ ^ Mqut- 

In the quark sector we choose a basis at the weak scale in which CKM rotation is 



entirely in the up-quark sector: we set the fermion rotation matrices ( A. 3D to be V 



Vuj^ = \^CKM and = ^d^i = 1- Note that this does not mean that remains 

diagonal at all scales - off-diagonal terms will be generated at higher scales due to RGE 
effects. Similarly, for charged leptons we set = Vg^, = 1 at Mz- 

Regarding the SUSY sector we work in the well-studied scenario specified by the pa- 
rameter set, 

mo, mi/2, ^0, tan/3, sgn{fj.) , (3.1) 



where GUT-scale boundary conditions are universal and defined by Eq. (|2.3| ). This choice 
of boundary conditions is frequently referred to as mSUGRA-like. Instead of scanning 
over the full parameter space which would be exceedingly computational intensive, we 
study specific points for each DM-allowed region. Throughout this work we take fi > 
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Figure 1: Code flowchart. Vi represent sparticles and Higgses arranged in the ascending order of 
their masses. The code is described in Appendix 



as suggested by measurements of the muon anomalous magnetic moment |56, ^] and 
set the pole mass of top quark mt = 171 GeV in accord with Tevatron data [^9|. For 
the DM relic density, we consider the conceptually simplest scenario in which the DM is 
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comprised only of the lightest neutralino Z\ that is thermally produced in the standard 
ACDM cosmology. We first calculate the neutralino Relic Density (RD) and LFV rates 
in our framework with SUSY-seesaw using points from Table || that have WMAP-allowed 
values in the mSUGRA framework without seesaw ~ a procedure commonly used in the 
literature. Then, if the RD turns out to be too high, as is frequently the case, we find new 
points consistent with the WMAP range ( [1.2| ) by adjusting SSB parameters and calculate 
the corresponding LFV rates. 

3.1 Large mixing 

We begin by considering the R^u = 1 case. Numerically we find that the GUT-scale 
Major ana mass matrix , 

Mat = diag (4.75 x 10"^ 4.75 x 10"^ l) • 1.398 x 10^^ GeV , (3.2) 

produces the spectrum mj^^ ~ 10^^ eV, rriiy^ — y/Am^ ~ 8x 10^^ eV and m^.^ ~ y^'Km^ ~ 
5 X 10~^ eV that is in accord with experimental limits (1.1). We chose the mass of the 



lightest RHN to be far above MsusY to prevent the unwanted mixing in the sneutrino 
sector. Equation ( |3.2| ) is in good agreement with our estimate ( |2.16D up to a factor of ~ 2 
reduction of up-quark Yukawa couplings (see for example Ref. [^]) due to the RGB effects. 
This is because fjy. Mat and the spectrum of light neutrinos are hierarchical and as such 
experience very little change in RGE evolution ||6l|| . 

In the Harrison-Scott parameterization ( 2.15| ), that we use for the mixing in the neu- 



trino Yukawa matrix at the GUT scale, the Dirac phase 6 and mixing angles are a function 
of (j) and X- S = -k/2 for </> = - see the left frame of Fig. |2|. It is known that LFV rates 
depend on the value of the unknown Dirac phase [^]. In addition, most LFV rates are 
very sensitive to the value of 6*13 [|l^] for which only an upper bound exists. To quantify 
these dependences, in the middle frame of Fig. ^ we show the branching ratio for /u ^ 67 
as a function of the Harrison-Scott parameters. We present them as enhancement factors 
relative to the rates at {(p, x} = {0, 0.294} for which 6*13 = 0.239 (or sin^ 6*13 = 0.056) and 
6 = 7r/2. We see that with ^13 fixed, the branching ratio changes by up to ~ 35% under 
variation of 5. The dependence on ^13 is greater and more complex: rates change by about 
two orders of magnitude for ^13 ranging from 0.239 to 0.045 (or sin^ 26*13 = 0.008), with the 
latter being the ultimate reach of the Daya Bay experiment |^3|. Closer to = x = 0, the 
rates change much faster and drop by several more orders of magnitude. The rates for the 
other 013-dependent LFV processes follow the same pattern as expected from Eqs. ( p^ ) 
and ( |2.9D . On the other hand, rates for r fij and r — >■ 3/U are relatively independent of 
^13 [0] • the right frame of Fig. ^ we show enhancement factor contours for r iJ.'j. We 
see that rates change with ^13 by only ~ 13%: this variation is an artifact of the paramter- 
ization ( 2.15| ), in which sin^ 612 = 0.33/ cos^ ^13. Also variation with respect to 6 is smaller 



- only up to ~ 10%. We numerically confirmed that r ^ 3/i rates behave similarly, as 
expected from Eq. (^). 

In Figs, m, ^ ^ and we show the LFV rates along with current experimental bounds 
(horizontal dashed lines) and projected future sensitivities (dash-dotted lines). To account 
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point A: large mixing, ff^^ =1, |i >0, =171 GeV 




X (rad) X (rad) % (rad) 



Figure 2: Dependence of the weak-scale Dirac phase (left), BR{ii ej) (middle) and 
BR{t — > fi-f) (right) on GUT-scale values of the Harrison-Scott parameters for benchmark point 
A with Yukawa unification parameter = 1. The contours for LFV branching ratios represent 
enhancement factors with respect to the point with S = tt/2, {(f), x} — {0, 0.294} that is marked 
by the blue cross; all angles are in radians. The thick red lines are iso-6'13 contours for ^13 = 0.239 
(at the CHOOZ bound ||6^; see Eq. LI) and for the ultimate 90% C. L. reach of the Daya Bay 
experiment ^13 = 0.045 | |63[ |. The contours remain essentially unchanged for all other points with 
large mixing and for — 3. 



for the aforementioned dependences on ^13 and 5, we show rates for {(j), x} = {0, 0.294} 
(for which sin^ ^13 = 0.056 and 5 = vr/2) as diagonally hatched bars and for (p = x = 
(resulting in sin^ 613 = 0) as solid bars. We also present rates for {(j), x} = {0, 0.022} 
(giving sin^ ^13 = 0.002 and 5 = 7r/2) as cross-hatched bars to indicate the upper limit on 



the rates if the Daya Bay |63| and Double Chooz [64| experiments produce a null result. 

In the bulk |2^ and the stau-coannihilation ||2^, ^ regions, the neutralino RD 
is within the WMAP range due to Zi interactions with the the lighter stau fi. Under 
universal boundary conditions fi is a dominantly right-handed state and as such remains 
unaffected by the neutrino Yukawa coupling^. Therefore, the use of mSUGRA values 
(points A and B) in models with RHNs still produces the correct RD. At point A in the 
R„u = 1 case, BR{fi 67) changes from 1.77 x lO"^'^ for = x = 0, to 8.57 x 10"^ 
for {(j), x} = {0, 0.294}; intermediate allowed (p and x values produce rates between these 
values. For r fij and r — > 3//, the dependence on (j) and x parameters is reversed: larger 
(f) and X produce smaller rates. For example, at point A in the Ri^u = 1 case, the t ^ fi^ 
branching fraction is 2.09 x 10"'^ for </> = x = and is 1.77 x 10"'' for {(j), x} = {0, 0.294}. 
We see that r fij and r ^ 3/i are excellent probes of LFV: the current experimental 



detailed discussion of neutrino Yukawa coupling effects on SUSY spectrum and DM observables can 
be found in Ref. Q. 
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Figure 3: Radiative LFV decay rates in the large mixing case for two values of R^u (c-f- Eqs- 2.4 



and 2.7). The heights of the solid bars show rates for exact tri-bimaximal mixing x — 4> — ^ 
at the GUT scale. Diagonally hatched bars represent {0, x} — {0, 0.294} yielding the maximum 
allowed ^13 and 5 = 7r/2. Cross-hatched bars represent {0, x} = {0, 0.022} which has 5 — 7r/2 
and 6*13 ~ 0.045 that is the ultimate reach of the Daya Bay experiment. Dashed lines represent the 
current bound and dash-dotted lines the projected sensitivity listed in Table ^. The letters denote 
various benchmark points presented in Tables H, ^ and ^. Points C through G have a RD above 
the WMAP bound. The corresponding WMAP-consistent points indicated by subscripts s, t and 
g are obtained by adjusting toqj ^0 and TOi/2j respectively. 

bound of r — > ^"f rules out the bulk region for all values of and x- 

In the stop-coannihilation region the picture is radically different. A naive use of 
input parameters for point C gives a neutralino relic density ^^^h? = 0.34, well above the 
WMAP bound. This is because ti is pushed to a higher mass causing stop-coannihilation to 
cease. To restore the stop-coannihilation mechanism and bring ^^^h? down to the WMAP 
range, one can counteract the effect of by adjusting SSB parameters. Adjusting the 
common scalar mass parameter mo, with the rest of SSB parameters held fixed, can lower 
the stop mass back to the desired mass range leading to a new WMAP-consistent point 
that we denote^ Cs with parameter values shown in Table ^. Increasing mo also makes 
sleptons lighter and increases their mixing, as can be seen from Eqs. (|2.4| ) and (2^), leading 



to a factor of ~ 2.8 increase in LFV rates as compared to point C. The stop mass can also 
be lowered by dialing the trilinear A-term, resulting in another modified point Ct- The 
rates increase with respect to values at point C by only about 10%. Alternatively, one can 
raise the Zi mass closer to that of the ti by adjusting the common gaugino mass parameter 
mx/2 producing a correct RD at the point Cg. Since the required increase is small, rates 



^Hereafter we use subscripts s, t and g to denote points obtained from those in Table |l| by adjusting the 
value of one the model parameters: scalar mass mo, frilinear coupling Aq or gaugino mass mi/21 respectively. 
Single and double primes are used to further distinguish modified points in the large mixing with R^u = 3 
and the small mixing with R^u = 1 scenarios. 
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Figure 4: Similar to Fig. || for t radiative LFV decays. For r fij, the {(p, x} = {0, 0.022} case 
is not shown as it produces rates very close to those for {(/>, x} — {0, 0.294} (cross-hatched). 



again grow only by ~ 10%. 



In the A- funnel region |24, 29], pushes the mass of the CP-odd Higgs boson A up 
and away from the resonance, resulting in a larger RD of 0.144 for point D. This can be 
reduced by decreasing the scalar mass parameter (point Dg) or increasing the value of the 
trilinear A-term (point Dt), which lower the A mass back to the resonance regime. In both 
cases LFV rates increase by about 8%. Increasing the gaugino mass parameter can increase 
the Zi mass to the resonance value = 0.5mA- But the A mass also grows with mi/2 
although slower than the neutralino mass, so a large dialing is required, resulting in point 
Dg. This large change in mi/2 increases the masses of the charginos in the loop resulting 
in about an order of magnitude drop in LFV rates. 

In the lower part of the HB /FP region [^] , at point E the RD is two orders of mag- 
nitude above the WMAP range. This is due to an increased value of // from the fj^ effect 
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Figure 5: Similar to Fig. |^ for trilepton LFV decays. 



that can be counteracted by increasing mo, resulting in new WMAP-allowed point Eg. The 
heavier sfermion spectrum causes LFV rates to decrease by about a factor of two. Adjust- 
ment of the trihnear parameter can somewhat lower /i Q , but not enough to get back into 
the HB/FP regime. It would be possible to reduce the RD by lowering mi/2 to 189 GeV, 
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Figure 6: Similar to Fig. | for LFV — > e conversion rates in titanium (upper frames) and 
aluminum (lower frames) targets. 

but at this value the chargino mass falls below the LEP2 bound of 103.5 GeV [|6^ ]. 

In the upper portion of the HB/FP region, the neutrino Yukawa couplings have an 
extremely large effect - the RD at point F becomes 12.3. Similarly to point E, the RD 
can be lowered by increasing the scalar mass parameter. In this part of the HB/FP region 
the chargino mass is sufficiently high that can be lowered without violating the LEP2 
chargino bound, resulting in a WMAP-allowed value at point Fg. Since charginos become 
lighter with this adjustment, LFV rates increase by about 40%. 

In the light Higgs resonance region [^, the neutrino Yukawa coupling also destabi- 
lizes the Zi annihilation mechanism producing too large a RD at point G. Although neither 
the Zi nor h masses are moved away from the desired regime 2m^^ = rrih, the resonance 
mechanism ceases to function because Zi becomes bino-like and can no longer couple to the 
higgs. The neutralino-higgs coupling can be restored by lowering n by increasing the scalar 
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Ds 
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51 
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51 
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i Z4 


n 
U 


oi 
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Es 
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10 
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Fs 
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10 
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Fg 
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10 
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Gs 


2435 
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-2000 


10 
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Gt 


2000 
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-1680 


10 
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Table 3: Modified benchmark points for mSUGRA-seesaw in case of the large mixing with R^u = 1 
obtained from their counterparts in Table |l| by adjusting the parameter highlighted in boldface to 
produce the RD dictated by WMAP. All dimensionful parameters are in GeV. For all points /i > 
and uit = 171 GeV. 



mass. This decreases LFV rates by about 20%. The desired value of the neutralino-higgs 
coupling can also be achieved by adjustment of the trilinear parameter resulting in point 
Gt- In this case, the LFV rates increase only marginally compared to point G. 



For R,yu = 3, we find that experimental limits (1.1) are satisfied for a GUT-scale 
Majorana mass matrix of the form, 

Mn = diag (4.5 x 10-^ 4.5 x 10-^ l) • 1.3 x 10^^ GeV, (3.3) 



which we use in subsequent computations. This is a simple rescaling of Eq. ( |3.2|) as expected 
from the seesaw formula. We have numerically verified that the dependence of the neutrino 
mixing parameters and LFV rates follow the same pattern shown in Fig. |2|. Thus, rates on 
Figs. are presented for the same choice of (j) and x discussed earlier. 

As for Ruu = 1, the RD in points A and B remain unaffected by the presence of 
additional neutrino Yukawa couplings. However, the very large neutrino Yukawa generates 
large off-diagonal terms in the slepton mass matrix, as can be seen from Eq. ( |2.4D , that 
boost LFV rates by more than an order of magnitude as compared to the R^u = 1 case. 
Nevertheless, this is still not enough to rule out point B for the whole range of x a-nd <j) 
values. 

At point C, the RD is too large: the ti mass is pushed further away from the Zi due 
to larger neutrino Yukawa effects. Thus restoration of the stop coannihilation mechanism 
requires larger adjustments of scalar mass and trilinear coupling parameters, leading to new 
WMAP-consistent points C'g and Cj listed in Table ^. Unlike the Ri^u = 1 case, adjustment 
of mi/2 cannot restore the stop-coannihilation: effects of f^, make ruf^ < m^^ and the RD 
is lowered to the WMAP range at mi^2 — '''25 GeV by the stau-coannihilation mechanism. 
Increasing mi/2 further makes fi the LSP before the stop-coannihilation regime can be 
reached. 
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Table 4: Similar to Table |^ for the large mixing and R^u = 3 case. 



In the A-funnel, larger neutrino Yukawas push uia away from the resonance resulting 
in ^^^h? = 0.33 at point D. Lowering the scalar mass parameter can bring the A mass back 
into the resonance regime at point D'^. The RD can also be lowered by either adjusting Aq 
to -692 GeV or raising mi/2 to 745 GeV. However, either of these bring the fi mass close 
to the Zi mass and activate the stau-coannihilation mechanism; further dialing of either 
parameter makes fi the LSP. 

In the HB/FP point E, fj, is pushed by neutrino Yukawa coupling to very large values 
resulting in large RD, ^^Jt? = 25. To compensate, one needs to dial the scalar mass 
parameter to very high values (point E'^). At such a large mo, sleptons become very heavy 
causing LFV rates to drop by about two orders of magnitude. Consequently, rates for LFV 
muon decay and fi — e conversion fall below current limits for all x (p values. For this 
point, fi is so large that the HB/FP regime cannot be recovered by adjusting or mi/2- 

Similarly, for point F we get an extremely high relic density ^^Ji^ = 72 that can 

be compensated by a very large scalar mass at point F'g. The LFV rates are lowered by 

a factor of ~ 15 so that muon LFV rates are below experimental bounds for all allowed 

mixing angles. As in the lower part of the HB/FP region, dialing or cannot bring 

0.7} h? back in accord with WMAP. 
^1 ^ 

In the /i-resonance point G, ^ is pushed so high that Zi becomes a pure bino state 
unable to couple to the Higgs boson. This can be compensated only by significantly in- 
creasing the scalar mass yielding a new WMAP-consistent point G'^ with LFV rates that 
are smaller by two orders of magnitude. 

3.2 Small mixing 



For the case of small mixing, we set iy{MQij'j') according to Eq. ( 2.17| ) and choose the 



neutrino spectrum to be = 6 x 10~^ eV, mj^j = 8 x 10^^ eV and iriy.j^ = 5 x 10~^ eV. Our 
choice for is constrained by the fact that we need Mtv^ ^ MsusY for our approximation 
to remain valid. In this scenario LFV rates do not depend on neutrino mixing angles: 
Eq. ( |2.17| ) fixes the neutrino Yukawa matrix completely and perturbations of the structure 
of Mjv do not produce significant changes in RGE evolution as can be seen from Eq. (p.^). 

For R^u = 1, fuh RGE evolution with our code yield the following eigenvalues of the 
Major ana mass matrix: 

Mn, ~ 8 X 10^ GeV , Mn, ^ 3.5 x 10^ GeV , MjVa ~ 2.5 x 10^^ GeV . (3.4) 
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This would appear to be in conflict with thermal leptogenesis, which requires the lightest 
Majorana mass to be heavier than about 10^ GeV |66]. Nevertherless, successful lepto- 



genesis is possible through the decay of the next-to-lightest RHN |67]. Also notice that 
M7V3 is closer to Mqut than it was in the case of large mixing. This, combined with small 
mixing in fj^, lead to rates that are several orders of magnitude smaller, putting them all 
significantly below current experimental bounds, as shown in Figs. ^, ^ and ^ 

As discussed in Section the use of mSUGRA values of points A and B still produces 
the correct RD. For example, at point A we have BR{fj, — > 67) = 1.1 x 10^^"^ that is barely 
above the reach of the future MEG experiment. In the ^-funnel, neutrino Yukawa couplings 
do affect the annihilation mechanism, but the effect is small and the RD remains within 
the WMAP range. For the other regions, adjustment of SSB parameters is necessary; the 
modified points are listed in Table ^. 

Even a decoupling scale as high 
as Mn^ ~ 10^^ GeV is enough to 
destabilize the stop-coannhilation 
mechanism, resulting in too large 
a RD for point C. Lowering the 
scalar mass parameter decreases the 
ii mass to the desired RD value 
at point Cg, leading to ~ 50% in- 
crease in the LFV rates. The stop 
coannihilation mechanism can also 
be restored by adjusting the trilin- 
ear A-term (point C") or the gaug- 
ino mass (point Cg). In both cases 

the LFV rates increase only slightly from those for point C. At all points /i — > 67 rates 
are slightly below the MEG reach, so this region will only be probed by // — > e conversion 
experiments. 
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Table 5: Similar to Tabic ^ for the small mixing and 
Ruu ~~ 1 case. 



At point E, the RD is also too high. For the reasons discussed in Section 3.1, only 
adjustment of mo is possible, leading to consistency with WMAP for the values at point 
Eg. Since the required shift is not as significant as in the large mixing case, the LFV rates 
drop only by ~ 30%. 

In the upper part of the HB/FP region at point D, the RD exceeds the WMAP value by 
two orders of magnitude. The desired higgsino content of Zi can be restored by adjusting 
the scalar or gaugino mass parameters, giving points and Fg . As a result, the LFV 
rates change by about 10% with respect to those at point D. 

At point F, the higgsino content of Zi is diminished by neutrino Yukawa RGE effects 
resulting in too high a RD. The /i-resonance mechanism can be restored by adjusting the 
scalar mass (point G'!.) or the trilinear A-term (point G'l). These adjustments change the 
LFV rates approximately by a factor of two with respect to the prediction for point F. 

For Ryu = 3, we find the eigenvalues of the Majorana mass matrix to be 



Mtvi ~ 7 X 10^ GeV , Mn^ ~ 3 x 10^° GeV , M^Vg ~ 2.2 x lO^*' GeV 



.16 



(3.5) 
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Figure 7: Radiative LFV decay rates in the small mixing case for benchmark points presented in 
Tables |l| and ^j. All rates are below current experimental bounds. Dash-dotted lines represent the 
projected future sensitivity listed in Table |^. 

The heaviest Majorana mass value is very close to Mgut ^ 2.3 x 10^*^ GeV. Because of 
this, the effect of on RD is neghgible and WMAP-consistent values are obtained for the 
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Figure 8: Similar to Fig. |^ for trilepton LFV decays. 

mSUGRA points in Table ||. Prom the resultant LFV rates in Figs. we see that the 
larger neutrino Yukawa coupling produces LFV rates that are smaller by almost an order 
of magnitude. This is opposite to what we saw in the large mixing scenarios where R^u = 3 
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Figure 9: Similar to Fig. |7| for LFV conversion rates in titanium (upper frames) and aluminum 
(lower frames) targets. 

rates were more than an order of magnitude greater that their Ryu = 1 cousins. This is 
also a direct consequence of the proximity of Mj^^ to Mqut'- the largest neutrino Yukawa 
decouples almost immediately and off-diagonal elements in the slepton doublet matrix are 
generated by the much smaller Yukawas of the first and second generations. 

4. Discussion 

We perfomed a detailed study of the LFV rates in the RD-allowed benchmark points, 
and demonstrated that the interconnection between the neutrino sector and neutralino 
dark matter is very important for predictions of LFV rates. Proper consideration of these 
effects change LFV rates by factor of a few to up to two orders of magnitude. We emphasize 
that although we used 50(10) models to set the structure of the neutrino Yukawa matrix 
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at the GUT scale, our results are generic; they hold in any type-I SUSY-seesaw scenario 
with large neutrino Yukawa couplings. 

The results in Section ^ imply the following about models with universal (or mSUGRA- 
like) SSB boundary conditions stipulated at the GUT scale: 

• The small mixing scenario is completely consistent with present experimental bounds 
on LFV. Upcoming fi ^ experiments will probe only a very small corner of param- 
eter space where both niQ and small and R^u = 1- Future ^ — > e conversion 
experiments, although suppressed by a factor ~ Za/ir with respect to /u — > 67, have 
better prospects due to the very well defined experimental signal. The PRIME ex- 
periment will be able to probe the entire bulk and stop-coannihilation regions as well 
as a significant portion of the ^-funnel region, while the Mu2e experiment will only 
be able to probe the bulk and stop-coannihilation regions for R^u = 1- 

• Contrary to naive expectations, in the small mixing case, order of magnitude smaller 
LFV rates are expected for Ri,^ = 3 than for R,^u = 1- Such small rates will only be 
probed by the PRIME /i ^ e conversion experiment with a Ti target. 

• Future — > 67 measurements will not rule out large mixing scenarios for any values 
of (j) and X in the mixing matrix of Eq. ( p.l5| ) because of the high sensitivity of this 
channel to ^13. On the other hand, if ^13 is close to the CHOOZ bound, only the 
HB/FP and /i-funnel regions are consistent with current limits. 

• The r ^ ;U7 channel is an excellent probe as it is not sensitive to ^13. Current 
experimental limits exclude the bulk region and the stop-coannihilation regions. For 
large Ruu, part of the A-funnel region is also excluded. Future experiments at Super 
Flavor factories |41| should be able to probe the A-funnel region almost entirely. 



• Trilepton decays are weaker probes due to a factor ~ a suppression of the rates as 
compared to the two body modes. Nevertheless, current data rule out ^13 close to 
the CHOOZ bound for some regions of SUSY parameter space. 

• /i — > e conversion in nuclei is the best probe of LFV. Future experiments will have 
sensitivity to almost the entire parameter space. In the large mixing case, /i — > e 
conversion is highly complementary to collider searches: it can probe large parts of 
the HB/FP region which can not be probed at the LHC. 

The upcoming Daya Bay and Double Chooz experiments will soon be able to probe 
013 independently of the Dirac phase down to sin^ ^13 = 0.002. A signal of nonzero ^13 will 
significantly reduce uncertainties in the predictions of LFV rates if the observed neutrino 
mixing arises dominantly from ii, as in the large mixing case. The 0i3-dependent LFV rates 
will be within about two orders of magnitude of the maximum values shown in Figs. |3[|^ 
thus further constraining model parameter space with current LFV data. For instance, the 
current ^ — > 67 bound rules out a significant portion of the A-funnel for the case of large 
mixing and either value of Ruu- With future — > 67 measurements we will be able to tell 
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if the type-I SUSY-seesaw can be realized in the stau-coannihilation region with large 
mixing regardless of R^u] see point B in Fig. ^. 

If the LHC finds a signal of SUSY, then GUTs become very appealing. One might be 
able to combine the knowledge of the sparticle spectrum from the LHC with results from 
LFV experiments to determine the value of R^u and/or get some information about the 
mixing pattern in the neutrino Yukawa. For example, if SUSY is found to be realized in 
the bulk region (point A), then a type-I SUSY-seesaw can only exist if the mixing in 
is small. Then measurements from the PRIME experiment could be used to identify Ruu- 
The situation becomes even more favorable if ^13 is known. For example, if SUSY is found 
to be consistent with the A-funnel region, then with the value of ^13 in hand, PRIME 
measurements will be able to test if the type-I SUSY-seesaw is operative for all mixing 
patterns and R,^u values. 



5. Conclusions 

In previous work [^] we demonstrated that large neutrino Yukawa couplings can signifi- 
cantly affect the neutralino relic density. This effect can be counteracted by the adjustment 
of SSB parameters, with concomitant changes in the low-energy phenomenology. In this 
work, we studied LFV processes in the type-I SUSY-seesaw properly taking into account 
neutrino Yukawa effects on the neutralino RD. For simplicity, we considered a scenario with 
flavor-blind universal (or mSUGRA-like) boundary conditions defined at the GUT scale. 
In the neutrino sector we utilized the "top-down" approach in which neutrino Yukawa and 
Majorana mass matrices are inputs at Mqut- We considered two cases for the neutrino- 



upquark unification parameter Ry^ (see Eq. 2.13 ) that are inspired by SO(IO) models. 
For each scenario we examined two extreme cases for the mixing in the neutrino Yukawa 
matrix. We found that the common practice of using WMAP allowed points of mSUGRA 
in models with RHNs overestimates the LFV rates in most regions of parameter space. 

In the Ruu = 1 case we found that the neutrino-neutralino interplay can result in 
significant changes in LFV predictions. The rates can change by a factor of up to 5 in 
the WMAP-allowed regions as compared to naive esimates. Effects are most prominent in 
regions with a large scalar mass parameter mo such as HP/FP and /i-funnel regions. If 
the mixing in the neutrino Yukawa matrix is small, then all LFV rates are below current 
experimental bounds. In the future, this case can be probed to some extent by the MEG 
experiment and by the PRIME and Mu2e conversion experiments. 

The case of very large unification parameter Ryu = 3, contrary to common lore, is not 
ruled out by current bounds on LFV processes even if the mixing in the neutrino Yukawa 
matrix is large. In the large mixing case a proper treatment of the neutralino-neutrino 
interplay leads to LFV rates that are smaller by about two orders of magnitude than 
naively expected. As result, many rates fall below current limits. Surprisingly, we found 
that if mixing in the neutrino Yukawa matrix is small, then for R^u = 3, the LFV rates are 
an order of magnitude smaller than for R^u = 1. If this scenario is realized, then only the 
future PRIME experiment will have sensitivity to some regions of the parameter space. 
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A. Notation and conventions 



Here we briefly list the most relevant equations to establish our formalism; we follow the 
notation and conventions of Ref. Q . 

The MSSM superpotential has the form 

/ = tJ-H^Hda + ^ ^{fu)ijeabQiH'^Uj + {id)ijQiHdaDj + {{e)ij^HdaEj , (A.l) 

where a, b are SU{2)l doublet indices, i,j are generation indices, eab is the totally anti- 
symmetric tensor with ei2 = 1 and the superscript c denotes charge conjugation. 
The soft SUSY breaking part of the Lagrangian is 



/-MSSM 
''-soft 



^ [MiAoAo + M2~XaXa + Ms^bOb] 



+ 



+ [bH^Hda + h.c] . 



(A.2) 



The gaugino fields {B = 1..8), \a = 1--3) and Aq transform according to the adjoint 
representations of 5C/(3)c, SU(2)l and U{1)y, respectively. 

We use the RL convention for the fermion mass term, Cmass = "(V'ij^-ijV'i + h.c), in 
which physical (real and diagonal) mass matrices read 



u V ur^m * Ul 



me = t^dVej^fJvt^ , (A.3) 

where Vu-, Vd are up- and down-higgs VEVs with v = \jv\ + v\ ~ 174 GeV. The unitary 
rotation matrices V, transform gauge eigenstates (unprimed) to mass eigenstates (primed) 
as follows: 



(VuL)iiULj , 
'^Li = (VdL)ii dLj , 



dm = (VdR)ijdijj , 
e'ffi = 0^en )ij eRj . 



(A.4) 
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The Cabibbo-Kobayashi-Maskawa (CKM) matrix |68| is ^ckm = ^ulV]j^. Due to differ 



ent matrix diagonalization conventions, our rotation matrices V, are hermitian conjugates 
of those in Ref. |9|,|7|. 



We work in the Super-CKM (SCKM) basis [71| where gluino vertices remain flavor 
diagonal. Here, the diagonalization of sfermion mass matrices proceeds in two steps. First, 
the squarks and sleptons are rotated "in parallel" to their fermionic superpartners 

d'Li = C^di^ )ij d-Lj , dffi = (VdH)ii dRj , 

1"^ = (y'eJijl'Lj , (A.5) 

where the SCKM scalar fields (primed) form supermultiplets with the corresponding fermion 
mass eigenstates, i.e., SCKM basis preserves the superfield structure after diagonalization 
of fermions. Next, 6x6 sfermion mass squared matrices in the SCKM basis are constructed: 

^ ( ^III + + D{uL)t -M?^^ + ^ cot /3 \ 
^ -MI{j^ + fi* cot P TUu M?^^ + ml + D{ur)1 J ' 



+ + -M|^^ + /.tan/?m/ 

^III + + D{eL)t -M?^^ + /itan/?me ^ 
-Mg[^ + fx* tan p rUe M?^^ + + D{dR)l^ 



(A.6) 



Ml = 

where D{f) are the hypercharge L)— term contributions to the corresponding sfermions, 1 



is the 3x3 unit matrix, mj are the diagonal fermion mass matrices of Eq. (A. 3), and the 
flavor-changing entries are contained in rotated SSB matrices, 

Mk = Vd,mi,vt^, M|^^ = Vd,mi,vt ^ , M|^^ = V^.a* , (A.7) 
MLl = Ve,m2 Vt, , M?^^ = Ve,m|vt^ , M?^^ = Ve,a*vt^ . 

Note that the squark doublet mass-squared SSB matrix hiq is rotated differently for Al| 
and A^j. Due to the absence of right-handed neutrino states in the MSSM, the sneutrino 
mass-squared matrix is a 3 x 3 matrix of the form. 

Ml = Ve^m^ Vt^ + D{9L)t . (A.8) 



Finally, the mass squared matrices ( [A.6| ) and ( [A.q ) are diagonalized to obtain sfermion 
mass eigenstates. These mass eigenstates are labelled in ascending mass order. 

To incorporate neutrino masses, we employ type-I seesaw mechanism where the MSSM 
is extended by three gauge singlet superfields N[ each of whose fermionic component is 
the left-handed anti-neutrino and scalar component is The extended superpotential 
has the form shown in Eq. (^]^). The neutrino Yukawa and Majorana mass matrices are 
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diagonalized in analogy with Eq. ( [A. 3D by unitary matrices Vj^j^ , V^^j^ and Vn according 
to 

rriD = Vu^^J^Wl^ , VnMn V?^ = diag(M^, , M^^ ,Mn,), (A.9) 

where m^) is the Dirac neutrino mass matrix. Apriori the eigenvalues Mjsi- are arbitrary, 
but the observed large difference between neutrino and charged lepton masses strongly 
suggests that the scale of Majorana masses M^jaj = max(M7V3) is much higher than the 
weak scale. Additional soft SUSY breaking terms should also be included so that the 
Lagrangian becomes, 



{a^)ijeabLiH''^i>l^j + ]r{h„)iji)Rii'jij + h.c. 



(A.IO) 



The eigenvalues of matrices and bjy and the square roots of the eigenvalues of m?^, are 
assumed to be of the order of the weak scale M^eak- 

Above the scale MMaj, the mass matrix for left-handed neutrinos is given by Eq. (p.2|). 
This is the minimal or type-I seesaw mechanism At low energies all RHNs decouple 
and the theory is governed by the effective superpotential, 

feff = Imssm + 2'^ij(^abLiHuedfLjHl . (A. 11) 

Here, k is a 3 x 3 complex symmetric coupling matrix that breaks lepton number explicitly, 
and is determined by matching conditions at RHN thresholds, 



(A.12) 



where (-^Nf, ) denotes the value as the scale of decoupling of the k-th generation RHN 
M^r^ is approached from above (below). 

After electroweak symmetry breaking, the Majorana mass matrix for light left-handed 
neutrinos is simply Adu = —kv'^. The symmetric 3x3 matrix Adu can be diagonalized as 

U^TW^U^ = , (A.13) 

where Uj^ is a unitary matrix and is a diagonal matrix of the physical neutrino masses 
mi, 7712 and ms. In labeling the (real non-negative) mass eigenstates we follow the usual 
convention that 1 and 2 denote states with the smallest mass-squared difference and mi < 
1712. The Maki-Nakagava-Sakata (MNS) mixing matrix of physical neutrinos [^] then has 
the form \mns = '^ej^'^t- Note that as a consequence of the seesaw mechanism, \J,y is in 
general different from the neutrino Yukawa diagonalization matrix Vt,j^ defined in Eq. ( 
The MNS matrix can be parametrized as 



V 



MNS 



( C12C13 S12C13 SlSe , ^ ^ 

-C23S12 - S23'Si3Ci2e"^ C23C12 - S23'Si3'Si2e"^ S23C13 | X diag(e*~ , 

\ S23S12 - C23'Sl3Cl2e*'^ -S23C12 - C23Sl3Sl2e*'^ C23C13 



(A.14) 
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with s 



sm 



9ij ctlld Cij 



COS 9ij, 9ij are mixing angles and the Dirac and Majorana 
CP-violating phases are 5, (/>i^2 S [0, 27r], respectively. 

The addition of the superfields also leads to an expansion of the sneutrino mass- 
squared matrix: it is now a 12 x 12 matrix, so that the relevant part of Lagrangian is 



C 3 nT 

2 



LtL 





(•^Lt 





2 

LtL 



-hi 



n . 



(A.15) 



where = (pj, u^^, P^, i^]^), is the 3x3 null matrix and 



(A.16) 



From this structure we see that there is sneutrino-antisneutrino mixing for right-handed 
states introduced by b^, with no corresponding terms in the left-handed sector. The Majo- 
rana mass matrix Mn contributes to the mass of the right-handed states and also results 
in the mixing of right-handed anti-sneutrino states with left-handed sneutrinos. Since Mn 
eigenvalues are much larger than the rest of the SSB parameters the matrix exhibits a 
seesaw type behavior, similar to the one for neutrinos: the 6x6 L-L block is of 0(-^^eaA:)' 
while R-L blocks are of 0{M^eakMMaj) and the R-R block is 0{Mli^-). Therefore the 
right-handed sneutrinos decouple and the phenomenologically relevant left-handed sneu- 
trinos have a mass-squared matrix of the familiar MSSM form (|A.8|). 



B. Description of ISAJET-M 

In this section we describe the 
algorithm used to perform the cal- 
culation. Our code, named ISAJET- 
M, is IS A JET 7.78 |5| modified to 
evolve all couplings and SSB pa- 
rameters, including the neutrino sec- 
tor, in full matrix form at the 2- 
loop level. The following Standard 
Model parameters are used as in- 
puts: fermion masses, the Z-boson 
pole mass M^, the fine structure 
constant a^^^{Mz), the strong cou- Table 6: SM input parameters [|| for ISAJET-M. All 
pling constant a^^{Mz) and CKM masses are in GeV and S^^^^ is in radians, 
angles in the "standard parametriza- 
tion" ; see Table ^. 



parameter 


value 


parameter 


value 


Mz 


91.1876 




171 


l/a*^^(Mz) 


127.918 


rrih - rUc 


3.42 


«F(Mz) 


0.1176 


nie X 10^ 


0.511 


sin^ 9wiMz) 


0.23122 




0.10566 


m„(2GeV) 


0.003 


nir 


1.77699 


md(2GeV) 


0.006 


sin Off ^ 


0.22715 


m,(2GeV) 


0.095 


sin^g^*^ 


0.04161 




1.25 


sin^fg^*^ 


0.003682 




4.20 


^CKM 


1.0884 
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In the first step a*^"^ and a^^'^ are evolved from Q = Mz down to Q = 2 GeV using 



.MS 



2-loop QCDxQED RGEs with the additional 3rd QCD loop |7|,|7|]. We use step-function 
decoupling of fermions at the scale of their running mass and include finite threshold 



corrections at 2 loops according to Ref. [76|. Then we compute the running lepton masses 



nil {2 GeV) from their pole masses mi using the 1-loop expression in the MS scheme 



a 



MS 



(Q) 



IT 



1 + - In ^ 



m. 



(B.l) 



Next, the two gauge couplings and all SM fermion masses (except for the top mass) are 
run up to Mz- Here fermion masses are converted from MS to DR scheme using formulae 
given in Ref. |7^, |78|: 



m^«(M^) = mf'{Mz) i^l 
rnf^iMz) = mf^{Mz) (l 



a 



29ai 



+ 



35i 



3tt 727r2 1287r2 19207r2 J 



a. 



29ai 



+ 



35i 



+ 



9l 



mf^(Mz) = mf^(Mz) ( 1 + 



MSt 



Sir 727r2 USn"^ 19207r2 J j^g 
35i 



(B.2) 



95L 

1287r2 6407r2 



MS 



For the lighter SM fermions these conversion corrections are neglected, i.e., we take m 



DR 

f 



m 



at Mz for f = u, d, s, e, /x. 
For the top quark, we obtain the running mass at Q = mt using the 2-loop QCD 



expression [79| 



mf^ijnt) = mt 



1 + 



5 as{mt) ( as{mt) 



+ 



vr 



TT 



-\2lo(yp 



(B.3) 



where Y^°°^ is the 2-loop piece. Note that we activate the top quark at the scale of its 
mass, so we have a 5-flavor scheme below Q = mt and all 6 flavors above it. 

to 



The obtained DR values of SM fermion masses are then substituted into Eq. (|^ 
calculate the Yukawa matrices at Mz in the gauge eigenbasis. A choice of basis is made by 



specifying fermion rotation matrices V, defined in Eq. ( A. 3 ). ISAJET-M also has options 
for performing calculations in the unmixed and dominant third- family approximations. 

The obtained weak scale values of gauge couplings and Yukawa matrices are evolved 
to the Grand Unification scale Mqut via MSSM RGEs in the DR scheme. The GUT 
scale Mgut is defined to be the scale at which gi = g2, where gi = ^hj^ig' is the hyper- 
charge coupling in the GUT-scale normalization. We do not impose an exact unification 
of the strong coupling (<^3 = g\ = 52) at Mqut, assuming that the resulting few percent 



discrepancy comes from GUT-scale threshold corrections |80|. 

For the evolution of gauge and Yukawa couplings, we use a multiscale effective theory 
approach proposed in Ref. , where heavy degrees of freedom are integrated out at each 
particle threshold. In Appendix we list the corrected Yukawa 1-loop RGEs from an 
unpublished erratum of Ref. . In the second-loop RGE terms we change from MSSM 



-27- 



formulae to SM expressions |6C, 7C] at a single scale, Q = MsjjSY = ^"^II "^Ir ' thus 
introducing an error of 3-loop order which is negligible. 

This "step beta-function approach" produces continuous matching conditions across 
thresholds. However, decoupling of a heavy particle also introduces finite shifts in RGE 
parameters a similar effect has long been known in QCD, where the decoupling of 
heavy quarks leads to shifts in the running masses of the light quarks |84|. Expressions 
for shifts induced by decoupling of each individual sparticle depend on the ordering of 
sparticle spectrum and are not yet known for the general case. Therefore, we implement 
these sparticle-induced finite shifts (to all three generations) collectively at a common scale 
Q = MsusY in the basis where Yukawa matrices are diagonal; we use 1-loop expressions of 
Ref. 1 85] without logarithmic terms that have already been resummed by the RGE evolu- 
tion. For the top Yukawa coupling additional 2-loop SUSY-QCD corrections are included 
according to Ref. [79|. These finite threshold corrections are particularly important for 
GUT theories since they change ratios of Yukawa couplings from those at the weak scale, 
as was emphasized in Ref. |86]. This multiscale approach is a generalization of the one 
used by the standard ISAJET, as is described in detail in Ref. and is preferred to 
single-scale decoupling when the sparticle mass hierarchy is large (as appears, for example, 
in the HB/FP region H of mSUGRA). 

At Mgut, the SSB boundary conditions are imposed and all SSB parameters along 
with gauge and Yukawa couplings are evolved back down to the weak scale Mz- For the 
SSB parameters we use 2-loop RGEs from Ref. [^ ] with the following conversion between 
notations: f, = Y,"^, a, = — h,"^, b = —B. Unlike the gauge and Yukawa couplings, 
where beta-functions change at every threshold, the SSB beta-functions remain those of 
the MSSM all the way down to Mz- We do not take into account threshold effects from 
the appearance of new couplings in the region of broken supersymmetry introduced in 
Ref. |7^. The entire parameter set is iteratively run between Mz and Mqut using full 
2-loop RGEs in matrix form until a stable solution is obtained. After each iteration the 
Higgs potential is minimized and sparticle/Higgs masses are recalculated. The obtained 
mass spectrum is used in the next iteration to appropriately account for sparticle threshold 
effects on the RGE evolution. 



The minimization of the RGE-improved potential is done using the tadpole method [B8|. 
We include 1-loop contributions from third generation sfermions, which dominate, as well 
as contributions from charginos, neutralinos and Higgs bosons. Computation is done at an 
optimized scale Q = MsusY, which effectively accounts for the leading 2-loop corrections. 
Since tadpoles that contribute to /i and b strongly depend on the parameters themselves, 
an iterative procedure is employed. Calculations of /i, b and tadpoles are iterated until 
consistent values with a precision of 0.1% are obtained; usually this requires 3-4 iterations. 

In the computation of sparticle masses we use SSB parameters extracted at their 
respective mass scales. Then, the SSB matrices are assembled and rotated to the SCKM 
basis as in Eq. ( A. 7 ). The resultant matrices are plugged into the sfermion mass-squared 



matrices ( A. 6 ). Instead of diagonalizing the full 6x6 matrices ( [A.6| ) we diagonalize three 
2x2 submatrices, thereby neglecting intragenerational mixings, which are required to be 
small by experimental limits on fiavor changing neutral current processes [Esf. For the 
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finite corrections, the full expressions of Ref. [^] for 1-loop self-energies are used. 

For the SUSY-seesaw we also use the multiscale approach that is mandatory to obtain 
correct values in the neutrino sector in the case of hierarchical RHNs . Our approach 



is identical to the one used in the REAP code |61|. Above the scale of the heaviest RHN 
we have a full MSSM+RHN setup, with the MSSM RGEs extended to include full 2- 
loop equations for f,^, a,y, and M^r |]l^, |l^]. As mentioned earlier, bilinear coupling 
matrix h,^ only introduces mixings for right-handed sneutrino states whose masses are 
0{MMaj) and is thus irrelevant for our analysis. We also take into account additional 
contributions to RGEs of ordinary MSSM parameters due to the RHN superfields up to 2- 



loop order [10, 15 , 61, For large neutrino Yukawa couplings these additional RGE terms 



cause changes to the MSSM sparticle spectrum that can have significant consequences for 



experimental rates |31, 32]. Below the scale of the lightest RHN, the RGEs are those of the 



MSSM, plus additional equations for the coupling k of the dimension-5 effective neutrino 



operator (A. 11) that are included at 2-loop level [Sl|. In the intermediate region, both k 
and some elements of and M^v are present. Transitions between them are done using 
matching conditions at RHN thresholds Eq. ( A.12| ). We match only at tree-level and neglect 



small finite threshold corrections. The position of RHN thresholds Mtv^. are determined 
by the eigenvalues of the RGE-evolving Major ana mass matrix M^v at those scales. Note 
that Eq. ( [A. 12] ) is valid only in the basis where M^v is diagonal at the threshold, which is 
different from the original basis at Mqjjt due to RGE effects. Below the scale Q = M^r^., 
the k—th. RHN superfields are absent from the theory so we remove the corresponding 
columns of and Sl^ matrices and the k—ih. row and column of m?^ and M^v in the basis 
where M^v is diagonal. 

Since some of the RGE parameters (gauge and Yukawa couplings) are defined at the 
weak scale, while the rest are set at Mgut, an iterative procedure is employed to solve the 
RGEs. Once a stable solution to the RGEs is obtained, the decay width and branching 
ratios for all the sparticles and Higgs bosons are calculated. This step as well as the previous 
computation of the mass spectrum is performed using standard ISAJET subroutines p5| . 
A graphical outline of our code, for a SUSY-seesaw model that uses top-down approach for 
the neutrino sector (i.e., f^, M^v are inputted at the GUT scale and physical light neutrino 
masses and mixing parameters are derived results), is shown in Fig. ^. 

ISAJET-M has facility to calculate and return the neutralino relic density Q.^J^^, 
neutralino-nucleon elastic cross sections, branching fractions for 6 — > S7 and Bs ^J'fJ'/T'^ 
decays, and supersymmetric contributions to the muon anomalous magnetic moment Aa^ = 



(g — 2)^/2. These computations are done, respectively, by the IsaReD ||2g], IsaReS |^ 



IsaBSG II, IsaBMM [|3| and IsaAMU |5| codes from the IsaTools package. In addi- 
tion, for scenarios with massive neutrinos, branching fractions for lepton flavor violating 
decays k — > Ijj and k 3lj as well as rates for — > e conversion in several nuclei can be 
calculated using expressions from Ref. |12|. 

C. Yukawa RGEs 

The following are the 1-loop RGEs for MSSM Yukawa coupling matrices with sparti- 
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cle/higgs thresholds from Ref. [Sl| with corrections from their unpubhshed erratum imple- 
mented: 



(4vr)2^i£k 
at 



^ ^ dt 



(JjLi 



k=l 



k=l 



k=l 



+- {c'e^ + s^Oh - ^c^Oh - Oh)) (fjf,*fj)., 

H€h [{-^^(^h + c') TV{3f:fJ} + c^iOh - 1) TV{3f,*fJ + CfJ}] 



+ 2 (c'^z. + ''^h) 



k=l 

N-, 



V d 

k=l k=l 

+(fj)., [.2(0, - 1) Tr{3f:fJ} + (c^e, + .2) TV{3qfJ + ^fj}] 



-if 



■d )ij 



^9t 



12 4 



36 



9^ + ^'^J?) 



l_ k=l k=l 

+(fj),, _ 1) Tr{3f:fJ} + (c20, + s^) TV{3f,*fJ + CfJ}] 

ror'ic;Q /I i \ ^ 



o 9 I 15 3 ^ 

T + 4^^ 



where s = sin a, c = cos a, a is Higgs mixing angle and the various 0-p's are equal to 
zero below the mass threshold of the respective particle and equal to one above it. The 
contributions from neutrino Yukawa couplings can be found in Refs. [10, 15]. 
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